DTRC/SHD- 1298-05  HOT-FILM  VELOCITY  MEASUREMENT  UNCERTAINTY  MARCH  1990 

FOR  DARPA  SUBCFF  EXPERIMENTS 


,3’ 

» - ^.vJL-X 

I  David  Taylor  Research  Center 

I  Bethesda,  MD  20084-5000 


AD-A220  002 

DTRC/SHD- 1298-05  MARCH  1990 
Ship  Hydromechanics  Department 

Technical  Report 

HOT-FILM  VELOCITY  MEASUREMENT  UNCERTAINTY 
FOR  DARPA  SUBOFF  EXPERIMENTS 

By 

James  N.  Blanton 
Thomas  J.  Forlini 
L.  Patrick  Purtell 


APPROVED  FOR  PUBLIC  RELEASE 
Distribution  Unlimited 


DTIC 

K  S3  ELEC!  e 
APR  021990 


MAJOR  DTRC  TECHNICAL  COMPONENTS 


CODE  Oil  DIRECTOR  OF  TECHNOLOGY,  PLANS  AND  ASSESSMENT 
12  SHIP  SYSTEMS  INTEGRATION  DEPARTMENT 

14  SHIP  ELECTROMAGNETIC  SIGNATURES  DEPARTMENT 

15  SHIP  HYDROMECHANICS  DEPARTMENT 

16  AVIATION  DEPARTMENT 

17  SHIP  STRUCTURES  AND  PROTECTION  DEPARTMENT 

18  COMPUTATION,  MATHEMATICS  &  LOGISTICS  DEPARTMENT 

19  SHIP  ACOUSTICS  DEPARTMENT 

27  PROPULSION  AND  AUXILIARY  SYSTEMS  DEPARTMENT 

28  SHIP  MATERIALS  ENGINEERING  DEPARTMENT 


DTRC  ISSUES  THREE  TYPES  OF  REPORTS: 

1 .  DTRC  reports,  a  formal  series,  contain  information  of  permanent  technical  value. 
They  carry  a  consecutive  numerical  identification  regardless  of  their  classification  or  the 
originating  department. 

2.  Departmental  reports,  a  semiformal  series,  contain  information  of  a  preliminary, 
temporary,  or  proprietary  nature  or  of  limited  interest  or  significance.  They  carry  a 
departmental  alphanumerical  identification. 

3.  Technical  memoranda,  an  informal  series,  contain  technical  documentation  of 
limited  use  and  interest.  They  are  primarily  working  papers  intended  for  internal  use.  They 
carry  an  identifying  number  which  indicates  their  type  and  the  numerical  code  of  the 
originating  department.  Any  distribution  outside  DTRC  must  be  approved  by  the  head  of 
the  originating  department  on  a  case-by-case  basis. 


NOW-OTNSHOC  5802/51  (R*v  2-M) 


■UNCLASSIFIED 


security  classification  of  this  paGe 


REPORT  DOCUMENTATION  PAGE 


la  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


2a  SECURITY  CLASSIFICATION  AUTHORITY 


2b  DECLASSIFICATION  /  DOWNGRADING  SCHEDULE 


4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

DTRC/SHD-1298-05 


lb  RESTRICTIVE  MARKINGS 


3  DISTRIBUTION /AVAILABILITY  OF  REPORT 
APPROVED  FOR  PUBLIC  RELEASE 
Distribution  Unlimited 


5  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


6a  NAME  OF  PERFORMING  ORGANIZATION  6b  OFFICE  SYMBOL  7a.  NAME  OF  MONITORING  ORGANIZATION 

David  Taylor  Research  Center  (lf  *pp,,m6/«) 

_  Code  1543 


6c  ADDRESS  (City.  Sutt.  s nd  ZIP  Cock)  I  7b  ADDRESS  (City.  Start.  and  ZIP  Code) 


8b  OFFICE  SYMBOL  I  9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  sppfiablt)  I 


Bethesda,  MD  20084-5000 


8a  NAME  OF  FUNDING /SPONSORING  1 8b  OFFICE  SYMBOL 

ORGANIZATION  I  (If  appliablt) 

Defense  Advanced  Research  Agenqy 


8c  ADDRESS  (City,  Stitt,  and  ZIP  Code) 

Submarine  Technology  Office 
STP  Support  Office 

1515  wilson  Blvd,  Suite  705,  VA  22209 


1 1  TITLE  (Incluck  Stcunty  Clitiificition) 

HOT-FILM  VELOCITY  MEASUREMENT  UNCERTAINTY  FOR  DARPA  SUBOFF  EXPERIMENTS 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM  I  PROJECT 
ELEMENT  NO  I  NO 


WORK  UNIT 
ACCESSION  NO 


12  PERSONAL  AUTHOR(S) 

James  N.  Blanton,  Thomas  J.  Forlini  and  L. Patrick  Purtell 


1 3a  type  of  report 

Departmental 


16  SUPPLEMENTARY  NOTATION 


13b  TIME  COVERED 
FROM _  TO 


14  DATE  OF  REPORT  (Year.  Month,  Diy)  IlS  PAGE  COUNT 

MARCH  1990  I 


COSATI  codes 


GROUP  SUBGROUP 


18  SUBJECT  TERMS  (Continue  on  rtvtrst  if  necessary  and  identity  by  block  numbtr) 


'9  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  numbtr) 


The  uncertainty  associated  with  measuring  the  three-component  velocity  for  the  DARPA  SUBOFF 
experiments  has  been  investigated.  Bias  uncertainty  due  to  temperature  variation,  probe  alignment,  analog 
to  digital  conversion,  speed  calibration  and  angle  calibration  resulted  in  an  overall  2#  bias  uncertainty  of 
±  221%  for  u/U^,  ±  2.28%  for  v/U,*,  and  £  1.79%  for  w/U„,.  Precision  (random)  2d  uncertainty 
depended  slightly  on  turbulence  intensity  levels.  For  low  turbulence  intensity  levels  (T1  «  0.5%), 
precision  uncertainty  was  ±0.04%  for  u/Uref,  ±  0.04%  for  vAJ^,  and  ±  0.06%  for  w/U*,.  For  high 
turbulence  intensity  levels  (TI  »  4%),  precision  uncertainty  was  ±10.28%  for  u/Uref;  ±  0.16%  for  v/U^, 
and  ±  020%  for  w/U^. 


20  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT  |2I  ABSTRACT  SECURITY  CLASSIFICATION 

□  UNCLASSIFIED/UNLIMITED  £3  SAME  AS  RPT  □  DTIC  USERS  I  UNCLASSIFIED _ 


22a  NAME  OF  RESPONSIBLE  INDIVIDUAL  |22b  TELEPHONE  (Include  A rtaCodt)  22c  OFFICE  SYMBOL 

James  n.  Blanton  I  (202)  227-1326  |  Code  1543 


DO  FORM  1473, 84  mar  *3  APR  edition  may  be  used  until  exhausted  SECURITY  CLASSIFICATION  OF  This 

All  other  edit, on,  are  obsolete  UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


CONTENTS 


NOTATION  .  v 

ABSTRACT .  1 

ADMINISTRATIVE  INFORMATION .  1 

INTRODUCTION .  1 

EXPERIMENTAL  FACILITIES .  2 

VELOCITY  MEASUREMENT  .  2 

HOT-FILM  ANALYSIS  .  2 

TEMPERATURE  COMPENSATION .  4 

PITOT-STATIC  VELOCITY  MEASUREMENT .  4 

BIAS  UNCERTAINTY  .  5 

TEMPERATURE  UNCERTAINTY  .  5 

PROBE  ALIGNMENT  UNCERTAINTY .  5 

PITOT-STATIC  VELOCITY  UNCERTAINTY .  7 

ANALOG  TO  DIGITAL  CONVERSION  UNCERTAINTY .  7 

SPEED  CALIBRATION  UNCERTAINTY .  7 

ANGLE  CALIBRATION  UNCERTAINTY  .  8 

SUMMARY  OF  BIAS  UNCERTAINTIES .  9 

PRECISION  UNCERTAINTY .  9 

SUMMARY  .  10 

ACKNOWLEDGMENTS .  11 

APPENDIX  A  THREE-COMPONENT  HOT-FILM  DATA  REDUCTION 

PROGRAMS  .  15 

APPENDIX  B  THREE-COMPONENT  HOT-FILM  SPEED  CALIBRATION 

PROGRAMS  .  23 

APPENDIX  C  THREE-COMPONENT  HOT-FILM  ANGLE  CALIBRATION 

PROGRAMS  .  31 

REFERENCES  .  39 


in 


FIGURES 


Page 

Figure  1.  Mean  and  turbulence  intensity  repeatability  measurements  (r/R=2.0) .  12 

Figure  2.  Mean  and  turbulence  intensity  repeatability  measurements  (r/R=0.25) .  13 


Figure  3.  Mean  velocity  repeatability  measurements  (r/R=0.72) .  14 

TABLES 

Table  1.  Typical  speed  calibration  data  and  errors .  8 

Table  2.  Bias  uncertainty  (%) .  9 


j  Accession  For  ; 

NTIS  GRA&I 

DTIC  TAB 

Unannounced 
Justification _ 

6 

□ 

By 

i 

Distribution/ 


|  Avail' 
! 

■Dlst 


mil  ity  Code 
’•nil  and/or 
Special 


m 


iv 


NOTATION 


■n 

E 

k 

L 

Qm 

AP 

r 

R 

T 

TI 

U,V,W 

UfOTAL 

v, 

x 

a 

♦ 

P 

e 

a 


direction  cosine  between  coordinate  direction  k  and  probe  direction  j 
bridge  voltage 
yaw  constant 
model  length 

effective  cooling  velocity  sensed  by  sensor  m 

pitot-static  dynamic  pressure 

probe  position  radius 

model  maximum  radius 

temperature  (°C) 

turbulence  intensity 

component  of  instantaneous  velocity  vector  in  direction  k  of  reference 
coordinate  system 

mean  velocity  in  three -component  directions  (x,  r,  0) 
total  velocity  vector  magnitude 

instantaneous  velocity  component  along  (parallel  to)  sensor  j 
axial  distance  measured  from  model  nose 
pitch  angle 

instantaneous  angle  between  total  velocity  vector  and  normal  to  sensor 
density 

angle  between  sensors,  or  probe  circumferential  position 
standard  deviation 


Subscripts 

s 

c 

e 

ref 


sensor 
calibration 
experiment 
reference  value 


v 


ABSTRACT 


The  uncertainty  associated  with  measuring  the  three-component 
velocity  for  the  DARPA  SUBOFF  experiments  has  been  investigated.  Bias 
uncertainty  due  to  temperature  variation,  probe  alignment,  analog  to  digital 
conversion,  speed  calibration  and  angle  calibration  resulted  in  an  overall  2 a 
bias  uncertainty  of  ±  2.21%  for  u/U^,  ±  2.28%  for  v and  ±  1.79%  for 
w/U„,.  Precision  (random)  2a  uncertainty  depended  slightly  on  turbulence 
intensity  levels.  For  low  turbulence  intensity  levels  (TI  »  0.5%),  precision 
uncertainty  was  ±  0.04%  for  u/Uref,  ±  0.04%  for  v/U^,  and  ±  0.06%  for 
w/Un,.  For  high  turbulence  intensity  levels  (TI  «*  4%),  precision  uncertainty 
was  ±  0.28%  for  u/Uref,  ±  0.16%  for  v/U^,  and  ±  0.20%  for  w/U^. 


ADMINISTRATIVE  INFORMATION 

The  work  described  in  this  report  was  funded  under  the  Defense  Advanced  Research 
Projects  Agency  (DARPA),  Task  Area  S1974-030,  Program  Element  63569N,  and 
performed  under  David  Taylor  Research  Center  (DTRC)  work  unit  1-1542-123. 

INTRODUCTION 

Experimental  velocity  measurements  were  performed  using  hot-film  anemometry  to 
assess  current  Computational  Fluid  Dynamics  (CFD)  capability  for  the  DARPA  SUBOFF 
Project  (Huang  et.  al,  1989').  Three  components  of  velocity  were  determined  from 
voltages  of  the  hot-film  sensors.  This  involved  analog  to  digital  data  conversion  and 
acquisition,  speed  and  angle  calibration,  temperature  compensation,  and  real-time  data 
reduction.  To  establish  a  level  of  confidence  in  these  velocity  measurements,  it  was 
necessary  to  estimate  the  experimental  uncertainty.  Bias  uncertainty  estimates  were  found 
from  perturbing  data  reduction  equations  with  estimates  of  individual  uncertainty  (e.g., 
temperature).  Precision  uncertainty  estimates  were  found  from  statistical  analysis  of 
repeated  measurements. 
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EXPERIMENTAL  FACILITIES 


The  David  Taylor  Research  Center  (DTRC)  Low  Turbulence  Wind  Tunnel  (LTWT) 
was  used  for  calibration  measurements  and  the  DTRC  Anechoic  Flow  Facility  (AFF)  was 
used  for  experiments  presented  in  this  report.  Velocity  measurements  were  normalized 
with  a  reference  velocity  measured  in  the  free-stream  at  an  axial  location  x/L  =  0.88, 
where  preliminary  measurements  showed  a  minimal  radial  velocity  gradient.  A  Datametrics 
Model  1174  Barocel  Electronic  Manometer  with  a  10  mm  Hg  pressure  transducer  was  used 
to  measure  the  dynamic  pressure  of  two  pitot-static  tubes  located  at  the  reference  plane. 

An  Omega  5830  thermistor  was  used  to  measure  free-stream  temperature.  A  TSI  IFA-100 
Intelligent  Flow  Analyzer  with  Model  150  anemometers  and  Model  157  signal  conditioning 
was  used  with  a  TSI  Model  12991-20  three  -component  hot-film  probe  to  measure  the 
instantaneous  velocity  components.  Data  from  the  Barocel  manometer  and  the 
anemometers  were  collected  using  a  Unix  based  Masscomp  5420  computer  with  16 
channels  of  333  khz  aggregate  analog  to  digital  conversion.  Compumotor  stepper  motors, 
controlled  over  the  IEEE-488  bus,  were  used  to  position  the  probe  during  both  calibration 
and  experimentation. 

Software,  written  in  the  ”C"  programing  language,  was  developed  on  the  Masscomp 
for  calibration,  data  acquisition,  probe  traversing,  data  reduction,  data  archiving,  and  initial 
plotting  of  velocity  (See  Appendices  A  through  C). 

VELOCITY  MEASUREMENT 

HOT-FILM  ANALYSIS 

The  total  velocity  vector,  Uto,^  can  be  related  to  the  effective  cooling  velocity, 

Qn>  (Hinze  (1975)J): 


2 


Qm  =  UNotal  (  cos2  <J>  +  k2  sin2  <j>  )  (1) 

For  three-component  analysis,  the  effective  cooling  velocity  as  a  function  of 
instantaneous  velocity  along  each  sensor,  Vj,  is  (Lakshminarayana,  1982s): 

3 

Q2m  =  2  Vs  (  kj  cos2  0mj  +  sin1  emj)  (2) 

where  k  is  the  yaw  constant. 

For  orthogonal  sensor  analysis,  0mj  =  90°  for  m  ^  j,  &  e-i  =  0°  for  m  =  j,  and  : 


k=l 


a2  =  1 


(3) 


where  a^H  are  the  direction  cosines  between  the  probe  coordinate  direction  and  the  sensor 
directions,  and: 

I  k2  1  1  | 

V2. 


Q2- 


k2  1  1 

1  k2  1 
1  1  k2 


(4) 


where  k=k,=kj=k3  is  assumed  for  identical  sensors. 

The  relationship  between  coordinate  system  velocity,  uk,  and  velocity  along  each 


sensor  is 


J 

Uk  =  2  a*,  Vj 

H 


(5) 


To  relate  the  anemometer  output,  E^,  to  the  response  of  the  sensor,  a  variation  of 
King’s  law  was  used: 

UJ*  =  A  +  B  E1,  (6) 

Separate  speed  and  angle  calibrations  were  performed.  For  speed  calibration,  U  = 
U^f  V  =  W  =  0  and: 


UJ*  =  A  +  B  E2 


(7) 
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A  series  of  voltages  for  various  reference  velocities  were  obtained.  An  iteration  scheme 
using  a  polynomial  least  square  fit  on  1/n  was  performed  to  obtain  n,  A,  and  B. 

For  angle  calibration,  the  probe  was  rotated  through  a  variety  of  yaw  and  pitch 
angles  in  the  uniform  free -stream  of  the  LTWT.  At  each  position,  three  components  of 
velocity  were  known  (determined  from  reference  velocity,  yaw  and  pitch  angles).  Using 
the  speed  calibration  found  in  (7),  and  substituting  into  (4),  the  cooling  velocity  was  found: 

Q2„  =  (A  +  B  E2r  [  1  +  (kJ  -  1)  aJ1>B  ]  (8) 

Once  Qm  was  calculated,  Vj  was  determined  from  (4),  and  then  uk  was  determined  from 
(5). 

From  a  set  of  possible  direction  cosines  (satisfying  restriction  of  (3)),  the  set  which 
results  in  the  lowest  difference  between  calculated  U  from  (5)  and  actual  velocity  from 
geometry  was  chosen  to  be  the  correct  direction  cosines. 

TEMPERATURE  COMPENSATION 


Since  the  heat  transfer  from  the  sensor  corresponds  to  velocity,  a  change  in  ambient 
temperature  results  in  a  change  in  measured  velocity.  This  can  be  compensated  for  if  the 
calibration  temperature,  Tc,  and  the  experiment  temperature,  Te,  are  known  (sensor 
temperature,  T#,  was  set  to  250  °C).  From  Buddhavarapu  (19864),  the  corrected 
anemometer  voltage,  Ec,  is  calculated  from  the  experiment  voltage,  Ee,  using: 


Ec2  =  E,:  x 


(  T,  -  Tc  ) 


(  T.  -  Te  ) 

Equation  (9)  can  then  be  substituted  into  (8). 
PITOT-STATIC  VELOCITY  MEASUREMENT 


(9) 


Reference  velocity  was  found  from: 


(10) 
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where  Ap  was  the  dynamic  pressure  from  the  reference  pitot  tubes.  Density  was  found 
from  a  fit  of  the  thermodynamic  properties  of  air  given  relative  humidity,  temperature,  and 
atmospheric  pressure  (Keenan  et.  al,  19805). 

BIAS  UNCERTAINTY 

The  bias  uncertainty  for  three-component  hot-film  measurements  can  be  divided  into 
five  categories: 

1.  Uncertainty  due  to  temperature. 

2.  Uncerta>  "  due  to  probe  alignment. 

3.  Uncertan^y  due  to  analog  to  digital  conversion. 

4.  Uncertainty  due  to  speed  calibration. 

5.  Uncertainty  due  to  angle  calibration. 

The  total  most  probable  bias  uncertainty  can  then  be  estimated  as  die  root  sum  square  of 
the  individual  uncertainties. 

TEMPERATURE  UNCERTAINTY 

The  accuracy  of  the  Omega  thermistor  is  stated  as  ±  0.3  °C  (±  0.2  °C  (instrument) 
+  ±  0.1  °C  (probe)).  The  uncertainty  in  velocity  obtained  by  perturbing  equations  (4),  (5) 
and  (8)  (using  the  method  of  Moffat,  1985‘)  by  ±  0.3  °C  was  ±  032%  for  uAJ^,  v/U 
and  w/U^. 

PROBE  ALIGNMENT  UNCERTAINTY 

The  three-component  probe  traversing  system  was  aligned  with  the  longitudinal  axis 
of  the  model.  Wake  surveys  were  performed  in  planes  normal  to  that  axis.  The  measured 
velocity  vectors  u,  v,  and  w  corresponded  to  components  in  a  model  coordinate  system,  x, 
r,  and  6.  Since  a  measurable  vertical  deflection  of  the  probe  tip  occurred  when  the 
traverse  was  moved  in  the  axial  direction,  alignment  was  performed  at  each  of  the  four 
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wake  survey  planes.  The  traversing  mechanism  was  adjustable  in  vertical  (pitch)  and 
horizontal  (yaw)  directions. 

For  the  two  survey  planes  on  the  body,  x/L  =  0.904  an  *  x/L  =  0.978,  the  traverse 
was  aligned  to  rotate  concentric  to  the  model  surface.  Pitch  of  the  traverse  was  adjusted 
with  measurements  of  the  distance  from  the  model  to  the  probe  at  0°  and  180°;  yaw  was 
adjusted  with  measurements  at  90°  and  270°.  These  measurements  were  performed  using  a 
scale  marked  in  1/64  in.  (0.04  cm)  increments.  The  probe  was  measured  concentric  to  the 
model  to  within  1/32  in.  (0.08  cm)  using  this  method;  alignment  is  estimated  to  be  within 
0.026°  of  the  model  longitudinal  axis. 

For  the  two  off-body  survey  planes,  x/L  =1.04  and  x/L  =  1.096,  the  traverse  was 
aligned  with  the  probe  stem  parallel  to  the  longitudinal  axis  of  the  model.  The  pitch  of 
the  traverse  was  adjusted  at  0°  using  a  line  level  on  the  probe  stem.  Previous 
measurements  using  a  small  laser  mounted  on  the  probe  stem  showed  this  method  accurate 
to  within  1/4  degree.  Yaw  of  the  traverse  was  adjusted  to  align  the  tube  parallel  to  the 
model  centerline  using  a  laser  aligned  along  the  pressure  taps  on  the  upper  surface  of  the 
model  and  a  vertical  rod  alternately  placed  on  the  forward  and  rearward  portions  of  the 
tube.  Alignment  within  0.3°  was  obtained  using  this  method. 

Based  on  a  nominal  reference  velocity  of  147  ft/sec  (44.8  m/s),  the  maximum  bias 
uncertainty  in  radial  and  circumferential  velocity  components  due  to  misalignment  in  pitch 
and  yaw  is  estimated  to  be  ±  0.5%  for  v/U^,  and  ±0.5%  for  w/U^.  The  bias  uncertainty 
in  axial  velocity  due  to  misalignment  is  negligible. 


6 


PITOT-STATIC  VELOCITY  UNCERTAINTY 


Measurement  of  reference  velocity  depends  on  temperature,  relative  humidity, 

atmospheric  pressure,  and  pitot-static  pressure  difference.  The  most  probable  uncertainties 

for  these  variables  are  estimated  as  follows: 

Temperature  ±  0.3  °C 

Relative  humidity  ±  5  % 

Atmospheric  pressure  ±  0.02  in  Hg 

Pitot-static  pressure  difference*  ±  0.01  mm  Hg 

The  uncertainty  in  reference  velocity  obtained  by  sequentially  perturbing  Equation  (10)  by 

each  of  these  values  and  then  root  sum  squaring  was  ±  0.1%. 

ANALOG  TO  DIGITAL  CONVERSION  UNCERTAINTY 

The  analog  to  digital  (A/D)  converter  resolution  was  12  bits  with  a  range  from  -10 

to  10  volts.  Therefore  one  A/D  count  was  20  volts  /  4096  counts  =  0.005  volts/count. 

The  anemometer  bridge  voltage  gain  was  6  and  the  A/D  gain  (on  the  A/D  board)  before 

digitizing  was  2.  Therefore,  to  determine  the  effect  of  one  A/D  count  on  velocity,  the 

bridge  voltage  was  perturbed  by  .005/12  =  0.0004  volts  (typical  reference  voltages  were 

used).  This  resulted  in  ±  0.07%  for  u/U^,  ±  0.13%  for  v/U^,  and  ±  0.12%  for  w/U^. 

SPEED  CALIBRATION  UNCERTAINTY 

Anemometer  voltages  for  various  reference  velocities  with  the  probe  normal  to  the 

free-stream  were  obtained  for  the  speed  calibration.  A  curve  fit  to  obtain  calibration  data 

reduction  equations  was  then  performed  (see  equation  (7)).  Calculated  values  of  U,  V,  and 

W  (as  distinct  from  the  true  reference  values,  U=U^,  V=0,  and  W=0)  were  found  by 

inserting  the  anemometer  voltages  into  the  calibrated  data  reduction  equations.  An  error, 


*For  Barocel  transducer:  ±  0.05%  of  reading  (8.6  mm  Hg)  ±  0.01%  range  (10  mm  Hg) 
±  1  A/D  count  (0.01  mm  Hg) 
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(Ur,)nllln1  -  U^/U^  was  obtained  for  each  calibration  velocity  as  shown  in  Table  1.  The 
uncertainty  (2a)  of  these  errors  for  a  typical  10  point  speed  calibration  was  found  to  be:  ± 
2.12%  for  u/U,^,  ±  1.76%  for  v/U,*,  and  ±  0.88%  for  w/U„,.  These  values  were  then  root 
sum  squared  with  pitot-tube  uncertainty  to  obtain  speed  calibration  uncertainty  of  ±  2.12% 
for  u/U,*,  ±  1.76%  for  v/U,*,  and  ±  0.88%  for  w/U„,  (pitot-tube  uncertainty  had 
negligible  effect). 


Table  1.  Typical  speed  calibration  data  and  errors. 


CALCULATED  REFERENCE  ERROR 


D 

V 

W 

(U 

-U~f>  /U„f 

V/U,., 

W/U,. 

ft /a 

ft/s 

ft/s 

1 

1 

ft/s 

1 

1 

% 

% 

% 

69.24 

-0.67 

-0.03 

1 

1 

68.95 

1 

1 

0.42 

O 

to 

00 

O 

O 

94.38 

-0.06 

0.03 

1 

96.41 

1 

-2.10 

-0.07 

0.03 

119.03 

0.56 

0.06 

1 

117.68 

1 

■ 

1.14 

0.48 

0.05 

125.94 

0.90 

0.22 

1 

124.17 

■ 

1.43 

0.72 

0.18 

130.66 

0.93 

0.41 

1 

129.89 

1 

| 

0.59 

0.72 

0.31 

136.84 

1.19 

0.29 

1 

136.35 

1 

0.36 

0.87 

0.21 

140.89 

1.34 

-0.49 

1 

141.54 

1 

t 

-0.46 

0.95 

-0.35 

145.82 

-0.81 

-1.57 

1 

146.99 

1 

| 

o 

GO 

O 

1 

-0.55 

-1.07 

152.86 

-1.08 

0.76 

1 

a 

152.41 

1 

| 

0.30 

-0.71 

0.50 

155.96 

-2.37 

0.34 

1 

I 

157.41 

1 

1 

-0.92 

-1.51 

0.22 

O 

1.06 

0.88 

0.44 

ANGLE  CALIBRATION  UNCERTAINTY 

Uncertainty  in  angle  calibration  can  be  calculated  similar  to  the  method  used  for  the 
speed  calibration.  However,  uncertainty  increases  with  increasing  flow  angle.  For  flow 
angles  25°  or  less,  2a  was  found  to  be  ±  2.08%  for  u/U,*,  ±  3.60%  for  v/U^,  arid  ±  3.80% 
for  w/U*,.  However,  flow  angles  greater  than  10°  were  not  observed  in  the  turbulence 
data.  Therefore,  2o  for  this  range  was  found  to  be  ±  0.52%  for  u/U^,  ±  1.30%  for  v/U^, 
and  ±  1.42%  for  w/U^. 
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The  uncertainty  in  determining  pitch  angle  a  in  the  angle  calibration  was  estimated 
at  ±  0.2°.  Perturbing  the  angle  calibration  data  by  0.2°  resulted  in  uncertainty  of  0.0%, 
0.01%,  and  0.2%  respectfully  for  u/U^,  v/U,*,  and  w/U^.  Root  sum  squaring  angle 
calibration,  a  and  pitot-tube  uncertainty  resulted  in  ±  0.53%  for  u/U„,,  ±  1.31%  for  v/U„,, 
and  ±  1.43%  for  w/U„,. 


SUMMARY  OF  BIAS  UNCERTAINTIES 
Root  sum  squaring  the  values  for  each  of  these  uncertainty  values,  as  presented  in 
Table  2,  resulted  in  an  overall  2 a  bias  uncertainty  of  ±  2.21%  for  u/U^,  ±  2.28%  for 
v/U^,  and  ±  1.79%  for  w/U„,  (Note:  the  uncertainty  due  to  A/D  conversion  is  negligible). 


Table  2.  Bias  uncertainty  (%). 


u/U„4 

v/UMf 

K/U„f 

« 

% 

« 

Temperature 

0.32 

0.32 

0.32 

Probe  alignment 

- 

0.5 

0.5 

A/D 

0.07 

0.13 

0.12 

Speed  calibration 

2.12 

1.76 

0.88 

Angle  calibration 

0.53 

1.31 

1.43 

Total  (root  a urn  square) 

2.21 

2.28 

1.79 

PRECISION  UNCERTAINTY 

Repeatability  measurements  were  made  to  determine  die  precision  (random) 
uncertainty  for  the  three-component  hot-film  measurements.  Measurements  were  made  both 
in  the  free-stream  away  from  the  model  wake  and  in  a  highly  turbulent  region  in  the  wake 
of  the  model.  A  single  measurement  at  each  position  required  <=  10-11  seconds  to  acquire 
data,  analyze  data,  and  move  to  the  next  position.  Typically,  1500  points  were  acquired  at 
150  Hz. 
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Figure  1  shows  mean  and  turbulence  intensity  versus  time  at  one  spatial  position 
(r/R  =  2.0,  x/L  =  1.096)  in  the  free-stream  (TI  =  0.5%)  for  thirty  points  taken 
consecutively.  Mean  and  2a  of  the  mean  velocity  were  0.955  ±  .04%  for  ii/Uref,  0.012  ± 
0.04%  for  v/U,*,  and  -0.011  ±  0.06%  for  w/U„,.  Figure  2  shows  profiles  for  one  spatial 
point  in  the  turbulent  wake  (TI  =  4%,  r/R  =  0.25,  x/L  =  1.04).  Mean  and  2a  of  the  mean 
velocity  were  0.653  ±  0.28%  for  u/Uref,  0.003  ±  0.16%  for  v/Ur<f,  and  -0.013  ±  0.20% 
for  w/Ur<r 

Figure  3  shows  mean  velocity  versus  time  at  r/R  =  0.72  and  x/L  =  1.04.  The  time 
involved  («  30  minutes)  is  representative  of  the  time  required  for  a  typical  constant  radius 
wake  survey.  In  the  time  between  8  and  25  minutes,  the  probe  was  moved  to  different 
radius  and  0  positions  and  then  returned  to  r/R  =  0.72  and  0  =  0°.  Mean  and  2a  of  the 
mean  (80  points)  velocity  were  0.938  ±  .16%  for  u/U^,  -0.008  ±  0.06%  for  v/U^,,  and  - 
0.020  ±  0.04%  for  w/U^,.  This  is  an  indication  of  the  repeatability  of  one  single  velocity 
measurement  over  the  typical  time  of  one  survey  and  the  repeatability  of  the  traverse  in 
returning  to  the  same  position. 


SUMMARY 

Bias  and  precision  uncertainty  associated  with  three  -component  hot-film  velocity 
measurement  have  been  investigated.  Bias  uncertainty  was  found  to  be  primarily  dependent 
upon  speed  calibration.  Other  significant  factors  were:  angle  calibration,  temperature,  and 
probe  alignment.  This  resulted  in  2a  bias  uncertainty  of  ±  2.21%  for  u/U„,,  ±  2.28%  for 
v/U„r,  and  ±  1.79%  for  w/U„,.  Precision  2a  uncertainty  was  found  to  have  a  slight 
dependency  upon  turbulence  intensity.  For  low  turbulence  intensity  levels  (TI  *»  0.5%) 
precision  uncertainty  was  ±  .04%  for  n/Uref,  ±  0.04%  for  v/U^,  and  ±  0.06%  for  w/U^. 
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For  high  turbulence  intensity  levels  (TI  «  4%),  precision  uncertainty  was  ±  0.28%  for 
u/Uref,  ±  0.16%  for  v/U,*,  and  ±  0.20%  for  w/U„,.  Therefore,  bias  uncertainty  was 
found  to  be  the  dominant  uncertainty. 
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Turbulence  Intensity  (X)  Mean  Velocity  /  Reference  Velocity 


Figure  2.  Mean  and  turbulence  intensity  repeatability  measurements  (r/R=0.25) 
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Figure  3.  Mean  velocity  repeatability  measurements  (r/R=0.72) 
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UILU 


APPENDIX  A 


THREE-COMPONENT  HOT-FILM  DATA  REDUCTION  PROGRAMS 

(Note:  These  programs  are  not  complete  in  themselves;  only  functions  (subroutines) 
pertinent  to  hot-film  data  reduction  are  included). 

/*  convert 3 . c  DTRC  Cod*  1543  7/89  */ 

/*  This  £11*  contains  th*  following  functions: 

volt3_to_y*loclty  () 
g*t_3c*l_d*t* ( ) 
allocate  () 
t*st_yolt  to  3v*l ( ) 

*/ 

(includ*  "hot3.h" 

PILE  *  apeed_cal ,  *dlr_c«l ; 

float  ‘slope, *incep, n, kappa; 
float  **dlr_cos, **inv_nat; 
float  cal_teap_c,  cal_teap; 

volt3_to_velocity (voltage, velocity, t*mp) 
doubl*  ‘voltage, ‘velocity; 
float  *t*ap; 

/*  Input:  voltag* [1,3] 

teaperatur* 

Return:  velocity [1, 3]  */ 

( 


convert  hf  voltag*  to  velocity 
gets  3  wir*  calibration  data 
Allocate  local  global  variables 
Test  voltag*  to  velocity  routine 


float  t*ap_eorr,  tensor,  sqr; 
float  *ps*odo,*g_aquar*d,  *v_v*locity; 
float  t*ap_c; 
int  1; 

/*  Allocate  sssory  and  unit  offset  variables  */ 
pseudo— vector (1,3) ; 
q_aquar*d— vector (1,3); 
v_veloclty— vector (1,3); 

/*  Convert  temperature  to  deg.  C  */ 
t*ap_o-  (  (*teep)  -32 . ) /I . 8; 

/*  Correct  voltages  for  temperature  */ 

/*  Calcinate  pseudo  velocity  (using  teaperatur*  compensation)  */ 
tesp_corr  -  eqrt ( (250 . -cal_t*ap_c) / (250 . -  teap_c) ) ; 

for(i-l;l<-3;l++)  ( 

sqr—voltag* [1] *t*ap_corr; 
teapor—  incep  [  1  ]  +  slope[l] *aqr*sqr; 
pseudo  [1]  -  pow  (teapor ,  n) ; 

) 

/•  Coapute  Q'a  given  aasuawd  direction  cosines  and  calibration  velocity  */ 
for (1*1; i<-3; 1++) 

q_squar*d [1 ] -pseudo [1 ] ‘pseudo [1 ] * 

(1.  +  (kappa* kappa  -  1.)* 
dlr_cos[l] [l]*dlr_cos(l] [1]); 

/*  Coapute  vj  by  sultipling  QA2  by  inverted  kappaA2  matrix  */ 
for(i«l;l<*3;i++)  ( 

v_velocity[i] -inv_**t [1] (i)*q_squsr*d[l]  + 
inv_mat[2] [1] *g_squar*d[2]  + 
lnv_aat  [3]  [1 j *g_aquar*d[3] ; 
if  (v_v*loclty [i]- <  0.)  ( 
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printf  ("**  Negative  V  **  "); 
v  velocity!!]-  -l.*v_velocity[i] ; 

)~ 

v  velocity [i)-aqrt (v  velocity [i] ) ; 

) 

/*  Compute  u,v,w  (confuted  velocity  in  probe  coordinates)  given  direction 
cosines  and  vj's  */ 

for(i-l;i<-3;i++) 

velocity tl]-dir_costl] [l]*v_veloeity[l]  + 

dir_cos[i] [ 2] *v_veloci ty [ 2 )  + 
dir_cos[ij [3] *v_velocity [3) ; 

velocity[2)  *-  -1.; 
velocity [3]  *-  -1.; 

/*  Deallocate  memory  */ 

£ree_vector (pseudo ,1,3); 
f ree_vector (q_aquared, 1,3); 
freevector (v_velocity, 1, 3)  ; 


/*  Function  to  allocate  memory  and  unit  offset  variables 
(Called  once  upon  entrance  to  main  program) 

NOTE:  These  variables  are  local  to  this  file  •/ 

allocate  () 

( 

slope— vector (1, 3) ; 
incep— vector (1,3); 
dir__coB-smtrix(l,3,l,3)  ; 
inv~mat -matrix (1, 3,1,3); 


/*  Read  in  direction  cosines  and  speed  calibration  data  */ 
get_3cal_dat a ( ) 

( 

int  i, j; 

float  a, in; 

char  flle_name[41] ; 

sprintf (file_naam, ’/home/blanton/hf /caldata/apdcal . %e",prabe_nua) ; 
speed_cal-fopen (flle_nssm,  "r") ; 
if  (lspeed_cal)  ( 

print f ("\n***  %a  could  not  be  opened\n", file_name) ; 
•*it(); 

) 

fscanf (speed_cal, "%f  %f ",  4n, 4cal_temp) ; 

print f ("\nnp  %6.4f  Calibration  tea p-  %4.1f",n,cal_teap) ; 
cal_te*p_o-  (cal_teag>-32 . )  /1 . 8 ; 

f or (i-1 ; i<— 3; i++)  ( 

fscanf ( speed_cal , "%f  %f", 4slope[i] , tineep[l] ) ; 

prlntf ("\nslope-  %6.4f  incep-  %6 . 4f" , slope [i] , incep [1] ) ; 

) 

printf ("\n") ; 
f close (epeed_cal) ; 

sprintf (filejnaaw, "/hoam/blanton/hf /caldata/dircos . %a", probe_num) ; 
dir  cal-f open (file  name,  "r") ; 
if  Ttdir_cal)  ( 

printf ("\n***  %s  could  not  be  opened\n”, file  name); 
exit(); 

) 

fscanf (dir  cal, "%f", 4 kappa) ; 
for (1-1; 1<— 3;i++) 

fscanf (dir_cal,"%f  %f  *f",4dir_eos[i) [1J , 

~  4 dir  cos [1]  [2] , 

6dir~cos[i] (3J); 


for  (i-l;i<-3;l++)  ( 
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for(j-l; j<-3; j++)  { 

print f  ("%6.4f  ",dir_co«[iJ  [  j] )  ; 

>  printf  ("\n") ; 

> 

fcloaa (dir_c»l) ; 

/*  Coaputa  Kffactivanaaa  Matrix  (Kappa  aquarad)  and  invart  */ 
kappa_matrlx (kappa, inv  Bat); 

) 

/*  Function  to  taat  voltaga_to_valocity  routina  */ 
taat  volt_to  3val() 

( 

doubla  voltaga [4] ,val[4] ; 
float  tarn p; 

printf ("\n\nlnput  3  valuaa  of  voltaga :  *); 
acanf(”%lf  %lf  %lf ", cvoltaga [1] , Cvoltaga[2] , tvoltaga[3] ) ; 
printf ("\Input  axpariaant  taaparatura:  ") ; 
acanf ("%f", stamp) ; 

printf ( "\nSl-%fi . 4f  E2-%6 . 4f  E3-46 . 4f " , voltaga [1] , voltaga [2] , voltaga [3] ) 
volt3_to_valodty  (voltaga,  val,  Stamp)  ; 

printf ("\n0  valocity  -  %f\nV  valocity  -  %f\nw  valocity  «  %f\n", 
val [1], val [2], val [3]); 

) 

/*  aatrix.c  DTRC  Coda  1543  7/89  */ 

/*  Thia  fila  containa  t 

kappa_matrix()  cosputa  affactivanaaa  (Kappa  A2)  matrix 
mat_inv()  invart  3X3  matrix  */ 

*induda  "hot3.h" 

/*  Cosputa  affactivnaaa  matrix  (kappa  **2)  */ 
kappajmatrix (kappa, inv) 
float  **inv, kappa; 

( 

float  **k_matrlx; 
int  i,j;  “ 
float  c2th,a2tb; 
char  tat [81]; 

k_matrlx  matrlx(l,3,l,3) ; 


for(i-l;i<-3;l++)  ( 

for (J-l; j<-3; j++)  { 

if  (i  I-  j)  c2th  -  0.; 
if  (i  I-  3)  a2th  -  1.; 
if  (i  —  3)  c2th  -  1.; 
if  (i  —  3)  a2th  -  0.; 

k  matrix [ij (3) -  kappa*kappa*c2th  +  a2th; 

) 

) 

/*  Invart  kappa A 2  matrix  */ 
mat_inv  (kjmatrix,  inv) ; 

/*  Fraa  mamory  */ 
fraa_matrix(k_matrix, 1, 3, 1, 3) ; 

» 

/*  To  Invart  a  3X3  matrix  */ 

/*  suit:  3X3  input  matrix 

inv:  3X3  Invart  ad  matrix 

dat:  dataradnant  */ 

mat_inv (mat, inv) 
float  **mat, **inv; 

( 
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float  dot; 
lnt  i,j; 

dot  -  B*t[l] [l]*mat[21  [21*mat[3J [3]  -  mat [1]  [1] *mat  [2J  [3J  *mat  [31  [2] 
-mat[l]  [2)*mat[2J  [l]*mat[3]  [3]  +  mat[ll [21*mat[21 [3]*mat[3] [1] 
+mat[lj  [3J*mat[2]  [U*mat[31  [2]  -  mat[U  [3]*mat  [2]  [2]*mat[3]  [11  ; 

inv[l][U  -  (mot [2]  [2J*mat[3]  [3]  -  mat  [2]  [31  *mat [3]  [21 )  /  dot; 
inv[2](2]  -  (mat  [11  [11 “mat [31  [3J  -  mat  [1]  [3]  *mat  [3][11 )  /  dot; 
inv[3]  [3J  -  (mattU[ll*mat[2][21  -  mat[ll  [2J*mat[2J  (11)  / 

/*  off  diagonal  olomonta  aro  transposed  wbon  calculated  */ 

lnv[ll [2]  -  -1* (mat [1] [2]  *mat [3] [3]  -  mat[U  [31*mat  [31  [21)  /  dot; 
lnv[U[3]  -  (mat[l]  [2]*mat[2]  [3]  -  mat  [1J  [3]  *mat  [2)  [2] )  /  dot; 

inv[2]  [1J  -  -l*(mat[2]  [ll«mat[3]  [3]  -  mat  [21  [3]  *mat  [3]  [1] )  /  dot; 
inv[2]  [3]  -  -1*  (mat  [1 J  [1]  *mat  [2]  [3]  -  mat[ll  [3]*mat[21  [1J )  /  dot; 

inv[3] [1]  -  (mat[2] [l]*mat[3] [2]  -  mat  [2]  [2J  *mat  [3]  [H )  /dot; 

inv[3] [2]  -  -1* (mat [1] [lj*mat[3] [2J  -  mat [11 [2] *mat [3) [1J )  /  dot; 


/*  otat.c  DTRC  Codo  1543  7/89  */ 

/*  This  fllo  contains  tho  following  functions: 

stats  ();  Conputo  moan  velocity,  rms  velocity, 

and  Reynold's  stress 

global  variables: 


Input: 

velocity  [num_data_j>ointa]  [3]  instantoous  u,v,w  velocity 
num_data-pointa  number  of  data  points 


Output: 

sman[l,3] 

rms[l,3] 

rstross(l,3] 

rstross[l] 
rstross [2] 
rstross[3] 

*/ 

# include  "hot3.h" 


moan  velocity 
rms  velocity 
Reynolds  stress 
(-u'v' )  (avg) 
(-u'w'l 
(-v'w'l 


stats  () 

( 

int  i, J; 

double  *sum,*sum2; 

/*  Allocate  memory  */ 
sum  rtvogtor(l,3) ; 
aum2— dvoctor  (1,3) ; 

for (i-l;i<«3; 1++) 
sum [ 1 1 “0 ; 

for(i-l;i<-num  data_points;  i++)  { 
for(J«i7 J<“3;3++)  ( 

sum[j]  +— velocity[iJ [ Jl ; 

1 

) 

for (i-1; i<“3; i++) 

mean [i J-sun[i] /num_data_polnts; 


for(i-l;i<-3;i++)  ( 

sum[l]»0; 

sum2[ll-0; 

for(l-l;K-f»um_data_points;i++)  ( 

for (J-l; J<“3;  J++)  ( 
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sum[j]  +-  (velocity [1] t j] -meant j] ) * (velocity!!) [j] -meant J) ) ; 

) 

■um2  [1]  +-  -1  .*  (velocity  [ij  [1]  -ra»an[l] )  *  (velocity (1)  [2]  -Man  [2] )  ; 

»um2 [2] +-  -1.* (velocity (i) [1] -maan[l] ) * (velocity [1] [3] -naan [3] ) ; 

»um2 [3  j  +“  -1 .  *  (velocityii  j  [2J  -waa[2] )  *  (velocity  [ij  [3] -mean [3) ) ; 

) 

for(i-l;i<-3;l++)  { 

ras[l)  -  sqrt (*um[i) /num_data_polnts) ; 
ratresa [i]-sum2 [i] /num_data_points; 

> 

/*  Deallocate  memory  */ 
f raa_dvact or ( aun, 1,3); 
f raa_dvactor ( sum2 ,1,3); 

) 

/*  aatup.c  DTRC  Coda  1543  7/89  */ 

/*  This  file  includaa  tba  following  function*: 

■*t_con*tant* ()  Road*  varloua  anvlornmental  constant*  from  flla 

"constants.dat",  prompt*  for  corractlon*,  aata 
currant  local  variables  and  writ**  corractlon* 
back  to  fila 

gat_dan*lty ( )  Raturn*  density  (slug*)  given  experiment  tamparatura, 

atmospharic  praaaura  and  daw  point  tamparatura 
af f_val 0  Return*  rafaranca  velocity  givan  pitot 

praaaura  difference  and  slug* 

ltwt_vel()  Return*  rafaranca  velocity  givan  pitot  praaaura 

difference, axparimant  tamparatura,  relative 
humidity  and  atmospharic  pressure 

*/ 

# include  "hot3.h" 

FXLB  * constants, **peed_cal; 

float  praaaure_atm; 
float  slugs, tempdew; 
float  raljhum; 

sat  constants  () 

{ 

char  temp[81] , test [81] ; 

/*  Input  constants  from  file  */ 

constant s-f open ( " /homa/blanton/hf /constants/constanta . dat " , "r+") ; 
if  (I constants)  { 

printf ("\n***  /hf /constants/constants. dat  could  not  b*  opened  ***\n") ; 

exit  () ; 

) 

facanf (constants, "%s", temp) ; 

if  (tolowar  (teap[0])  —  • x ')  proba_typa  -  XWIRX; 
if (tolowar (taap[0j 1  —  *t')  probe_type  -  TRIWIRB; 
if  (tolowar  (temp  [0J>  I-  'x'  cs  tolowar  (tamp  [0]  >  I-  't')  { 
print f("\n  Probe  type  varlabla  not  correct \n") ; 
exit  () ; 

) 

facanf (constants, "%a",prob*_nun) ; 
facanf (constants, "%a", tamp) ; 

if  (tolowar  (teop[0J)  —  'a')  wind_tunnal  -  AFT; 
if (tolowar (tamp [0])  —  '1')  wind“tunnal  -  LTWT; 
if  (tolowar  (taop[0])  t-  '1'  cc  tolowar  (taap[0) )  !-  'a')  ( 
printf ("\n  Wind  tunnel  variable  not  correct\n") ; 
axlt(); 

) 

facanf (constants, "%f " , spreaaure_atm) ; 
if  (wind_tuimel  —  AFP  )  { 
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print f ("\n\t\tYou  art  using  AFF  wind  tunnsl  option\n")  ; 
f  scanf  (constants,  "%f  ",  4t8iup_dsw)  ; 

) 

if  (wind_tunnsl  —  LTWT)  { 

print f ("\n\t\tYou  ars  using  LTWT  wind  tunnsl  optlon\n"); 
f scanf (constants, "%f ", Srsl_hujn) ; 

) 

do  ( 


/*  print  constants  to  scrssn  */ 
if (prob*_typ*  — -  XWIRE  ) 

print f ("\n\nl .  Probs  Typs  ■  XWIRE") ; 
if (prob*_typ*  —  TRIWIRB  ) 

printf ("\n\nl.  Probs  Typs  -  TRIWIRB") ; 
printf  ("\n2.  Probs  Numbs r  -  %s",probs_num)  ; 

printf ("\n3 .  Atmosphsric  prsssurs  -  %5.2f  in.  Hg",prsasurs_atm) ; 
if  (wind_tunnsl  —  AFP  ) 

printf ("\n4.  Dsw  Point  Tsmp  -  %4.1f  dsg  F",  tsnp_dsw) ; 

if  (wind  tunnsl  — 1  LTWT) 

printf ("\n4.  Rslativs  Humidity  -  %4.1f%%",rsl_hum) ; 
printf ("\nq.  Contlnus  and  Savs\n\n:  "); 

gsts(tsat) ; 
switch  (tsst[OJ)  ( 

/*  Chang*  constants  (if  nscsssary)  */ 
csss  '1' : 
do  { 

printf (*\nKntsr  Probs  Typs\n\tx  -  X  Film\n\tt  -  3  Conponsnt\n:  ") 

scanf ("%a", tsmp) ;  _ 

if (tolowsr(tsop[OJ)  —  'x')  probs_typs  -  XWIRE; 
if (tolowsr(tsmp[OJ )  —  't')  probstyps  -  TRIWIRB; 

)  whils  (tolowsr  (tsmp  [  0  ] )  I-  'x'  fit  tolowsr  (tsmp[0) )  I-  't'); 
gstchar ( ) ;  /*  Ost  carriags  rsturn  lsft  ovsr  from  scanf?  */ 

cass  '2' ! 

printf <"\n\n\tl96  -  Thrss  Conponsnt  Sarial  *#1196*); 
printf ("\n\tt89  -  Thrss  Cospnsnt  Ssrial  #91089"); 
printf ("\n\ta30  -  Two  Cosponsnt  Ssrial  ♦???"); 
printf ("\n\tx51  -  Two  Conponsnt  ssrial  #?»?"); 
printf (*\nBntsr  Probs  Numbsr:  "); 
scanf ("%s",probs_num) ; 

gstchar ();  /*  Ost  carriags  rsturn  lsft  ovsr  from  scanf?  */ 

brsak; 
cass  '3' t 

printf ("\nBntsr  Atmosphsric  Prsssurs  (in.  Bg) :  "); 
scanf ("%f ", «prsasurs_atm) ; 

gstchar ();  /*  Ost  carriags  rsturn  lsft  ovsr  from  scanf?  */ 

brsak; 

cass  ' 4' > 

if  (wind_tunnsl  «*■  AFF)  { 

printf  ("\nBntsr  Dsw  Point  Tsopsraturs  (dsg  F)  >  ”); 
scanf ("%f ", ttsmp_dsw) ; 

) 

if  (wind  tunnsl  —  LTWT)  ( 

printf ("\n*ntsr  Rslativs  Humidity  (**) i  "); 
scanf (■%£", firaljhum) ; 

gstchar ();  /*  Ost  carriags  rstum  lsft  ovsr  from  scanf?  */ 

brsak; 
cass  '  q'  i 

/*  savs  constants  to  fils  */ 

rswlnd (constants) ; 

if (  a  bs_typs  — -  XWIRE  ) 

fprintf ( constant s,"XWXRB\n") ; 
if (probs  typs  —  TRIWIRB  ) 

fprintf (constants , "TRIWlRB\n" ) ; 
fprintf (constants, "*s\n",probs_num) ; 
if  (wind_tunnsl  —  AFF)  ( 

fprintf (constants, "AFF\n") ; 

fprintf (constants, "%3 . 2f \n" , prsssurs_atm) ; 
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fprintf (constant*, "%4 . If \n", tomp_dew) ; 

) 

if  (wind_tunnol  —  LTWT)  ( 

fprintf  (constant*,  "LTWT\n") ; 

fprintf (constants, "%5 . 2f\n”,pre**ure_atm) ; 

fprintf (constants, "%4 . If \n" , rol_hum)7 

) 

fdoso  (constants) ; 
break; 

) 

)  while  (tost [0]  I-  'q'); 


) 

/*  CALCULATE  DENSITY  */ 

/*  Rotum  dansity  — >  slugs 

taag>  (temperature)  — >  dog  F  * / 

float  get_den*ity (tamp) 
float  tonp; 

( 

float  acrl,*cr2,  slugs; 

/*  Curve  fit  from  Dave  Fry  */ 

/*  Input  dog  F  */ 

serl  «  .331 0  -  0.01256  *  tomp_dew  +  2 . 573e-4*teap_dew*teap_dew 
scr2  “  491.63/(459.63  4  tomp)  *  (2 . 54*proasuro  atm  -  scrl)/76. 
slug*  -  *cr2  *  1.2929o-3  *  1.9403; 
rotum  (slugs) ; 

) 

/*  CALCULATE  REFERENCE  VELOCITY  FOR  AFF  KIND  TUNNEL  CONDITION  */ 

/*  Rotum* :  velocity  ->  ft/s 

pitot_pross  difforonco  ->  sn  Eg 
taf  (temperature)  ->  dog  F  */ 

float  aff_vol (pitot_pro*s, ta) 
double  pitot_proas; 
float  ta; 

< 

float  rho,rofvol; 

rho— got_don*lty (ta) ; 

refvel-aqrt (5 . 569*pitot  jproas/rho) ; 

rotum  (rofvol) ; 

) 

/•  CALCULATE  REFERENCE  VELOCITY  FOR  LTWT  WIND  TUNNEL  CONDITION  */ 

/*  Returns:  velocity  •>  ft/s 

pitotjpross  difforonco  ->  mm  Hg 
taf  (teaperature)  ->  dog  F  */ 

float  ltvt_vol (pitot_pross, taf) 
double  pitot_press; 
float  taf; 

( 

float  ta; 

float  ps,pv,pa,w,  r,  gassoa,  tor; 
float  refvel(pa_am_hg; 

pa_mm_hyproa*uro_atm*25 . 4 ;  /*  convert  in.  Bg  to  ant  Hg  */ 

/*  Curve  fit  and  calculation  from  F.  Furtell  */ 
ta- (taf-32. ) /I . 8; 

ps— 3 . 049- . 08036*taf44 . 802e-3*taf *taf-3 . 536o-5*taf*taf»taf 
+4 . 14e-7*taf *taf *taf *taf ; 

pv-rel_hu**ps/100 .  ; 

pa-pa_ao_hg-pv ; 

w-0 . 62188*pv/pa; 

r»(w*2759.  92  +  1716.35) / (1 .4w) ; 

ga asaa-  (0 . 239840 . 4446*w)  /  (0 . 171240 . 3343*w) ; 

tor—(ta4273. 18) *1 . 8; 

rofvol-sqrt  (2 .  *r*tor*pltot_pre»*/pa_smi_bg)  * 
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) 


(1 .  +pitot_prea*/pa_imiiJhg/4 .  /guana.)  ; 
return (refvel) ; 


/*  hot3.h  */ 

♦include  <atdio.h> 

♦  Include  "/dtrc/indude/nrutll  .h" 
♦Include  <math.h> 

♦define  PI  3.1415926 


/*  Local  parameter*  (note  explicit  typing) :  */ 

♦define  NOLL  0 
♦define  SIZBOFOCA  1000 
♦define  NOLLINTRAY  -2 

char  probe_num[4] ;  /*  Number  for  apecific  probe  */ 

float  ** velocity;  /*  instant  velocity  array  */ 

short  *raw_lnt_voltage;  /*  Raw  data  array  */ 

float  r*f_veloclty;  /*  free  stream  velocity  »/ 

float  *asan, *ras, *rstre*a;  /*  3  deep  vectors  for  each  ehn.*/ 

/*  float  •ean[4),rsw[4],rstreas(4];  */ 

int  num_data_pointa;  /*  Number  of  points /channel  to  be  taken  */ 
enum  (  AFF  ,  LINT  )  wind_tunn*l; 
enum  {  XNIRS  ,  TRZN1RB  }  probe_type; 

int  pathno;  7*  path  number  for  IKK  406  */ 

float  n»an_vel[300] [3] , rm*_vel [300] [3] , re_streas [300] [3]; 
float  plot_rad[300] ,  plot_theta[300] ; 
int  group_num_array[300] ; 
int  theta_addreaa,rad_address; 

/*  struct  shared  { 

int  num_of_plot_pta; 
char  cosuentl  [89] ; 
char  consent 2  [89] ; 
float  plot_rad[300] ; 
float  plot_theta[300] ; 
float  aaan_vel(300] [3] ; 
float  ih  vel  [300]  [3]; 
float  re  stress [300] [3] ; 

In¬ 
struct  shared  *data;  */ 
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APPENDIX  B 


THREE-COMPONENT  HOT-FILM  SPEED  CALIBRATION  PROGRAMS 

(Note:  These  programs  are  not  complete  in  themselves;  only  functions  (subroutines) 
pertinent  to  hot-film  speed  calibration  are  included). 

/*  ap**dcal.c  DTRC  Cod*  1543  5/89  */ 

/*  THREE  COMPONENT  HOFTFILM  CALIBRATION  MODULE 

Tbl*  program  lnputa  3  hot  £11*  voltage*  and  1  pitot  tub*  voltag* 

for  a  rang*  of  tunn*l  v*lociti*a  and  compute*  alop*a  and  int*rc*pta  for  *acb  velocity  consonant 
*/ 

♦lnclud*  "cal .h" 

float  cal_t*ng>[30] ; 
char  f ll*_nam* [ 41] ; 

PILE  *cal_fil*; 

doubl*  hf_voltag*_m*an[30] [3] ; 
doubl*  pitot_v*locity_»*an[30] ; 
float  avg_eal_t*mp; 

maln() 

I 

lnt  don*,  1,  j,  cat; 

char  cb«ck[81] , lchk[81] , t*at*r[81] ; 

doubl*  h£_avg_volt [30] [3] ;  /*  dlff*r*nt  from  angvel.c  */ 

float  pltot_avg_v*l [30) ; 
float  alop*[3] , inc*p[3 j ; 

float  g*t_t«np();  /*  function  declaration*  */ 

float  ltwt_v*l ( ) , af f_v*l ( ) ; 

float  p_volt ; 

double  aum[4J; 

float  t*ag>_aax; 

float  t*ap_mln; 


/*  Allocate  memory  and  unit  offaet  variable*  */ 
hf_volt-dmatrlx  (1,  nframaa,  1, 3) ; 
pltot_volt«dv*ctor (l,nfram*a) > 

/*  8*t  *nvl orman t  */ 

a*t_conatanta  ()  ; 

don*  ■  0; 

do  ( 

prlntf ("\n\n\nlnt*r  :\n\ 

\tl.  Take  Calibration  Data\n\ 

\t2 .  Analya*  Calibration  Data\n\ 

\t3.  Chang*  Conatanta\n\ 

\tq  Quit\n\n\ 

\t : ; 


acanf("%e",  chock); 

awltch  (tolow*r (check [0] ) )  { 

caa*  '1'  : 

/*  S*t  •nviorment  */ 
a* t  atod  conatanta  () ; 

/*  Initialia*  US  488  */ 
inlt_gpib()  ; 

/*  Open  output  fil*  */ 

aprlntf (fll*_nam», "/hom*/blanton/hf/c*ldata/raw*p**d. %a",prob*_num) ; 
cal  fil*-fop*n(fil*  name,  "r"); 

“if  <lcal_file)  ”< 
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printf ("\n***  Will  craata  rawapaad . %a\n" , proba  nun); 

> 

if  <cal_fila)  ( 

printf("\nDo  you  want  to  appand  to  rawapaad.%a?  ",proba_num) ; 
acanf ("%a", taatar) ; 
if (tolowar (taatar [0] )  I-  'y'(  ( 

fcloaa (cal_fila) ; 
fopan(fila  nama, "w") ; 

} 

) 

fcloaa (cal_fila) ; 

tamp_max— 0; 
tamp~min-100 . ; 
ent-O; 


do  { 

printf("\n  Taking  Data  Run  #%d\n\n",cnt  +  1) ; 

gat_ang_voltag# ( )  ; 
cal_tamp[cnt ]  —  gat_tamp  () ; 

printf  ("\n  Calibration  Taniparatura  ia  %f\n", cal_tanp[ent] ) ; 

if  (cal_tamp[cnt]  <  tamp_m±n)  tamp_min— cal_tamp"fcnt]  ; 

if  (cal_taBp[cntj  >  tamp_max)  tanp_max-cal_tamp[cnt] ; 

for  (  i-0;  i<-2;  i++) 

aum[i]  -  0.; 

for  (  i-0;  i<— 2;  i++>  { 

for  (J-1;  j<-nframaa;  J++)  ( 

aum[i J  +-  hf  volttJJ  [i+1]  ; 

) 


for  (  i-0;  i<-2;  i++)  ( 

bf_avg_voltfcnt] fi]  -  aum(i]  /  nframaa; 
print fT"Voltaga  #%d  ia  %f\n",i, aum[i] /nframaa) ; 

) 

•um[0]  -  0.; 

for  (J-1;  J<-nfra*aa;  J++  )  ( 

aua[0)  +-  pitot  wolt[JJ; 

} 

p_volt  -  aua[0]  /  nfrUMi) 


if (wind_tunnal  —  LTWT) 

pitot_avg_val tent] -ltwt_val (p_yolt, cal_tanp [ent ] ) ; 
if (wind_tunnai  —  ATT) 

pitot_avg_ral [ent]-af f_val (p_volt , cal_tanp[cnt ] ) ; 


aprlntf (fila_nama, "/boaM/blanton/bf/caldata/rawapaad.*a",proba_nun) ; 

cal  fila-fopan(fila  nna,  "a"); 

f printf (cal_fila, "%Io . 3f  *10 . 3f %10 . 5f %10 . 4f *10 . 3f \n" , 
hf_avg_volt tent] (0J , 
hf~arg~Tolt [ent ] tlj, 
hf_awg_volt(cnt] [2] , 
pitot_avg_ral tent] , 
cal_taap[cnt] ) ; 
fcloaa (caljfila) ; 

printf ("\nValocity  ia  %f\n",pitot_avg_val tent] ) ; 

printf ("\n\nHit  c  to  taka  naxt  valocity  or  'q'  to  continual  ") ; 

acanf ("*a", lchk) ; 

cnt++; 

)  wbila  (tolowar (lchk [0])  I-  'q'); 
auai[0]-0; 

for (i-1 ; i<cnt ; i++) 

aun [  0  ]  +-  cal_taap[i]  ; 
avg_cal_taop— aum [ 0] / (cnt-1) ; 

printf ("\nkvg  tamp-%4 . If  Min  Taa^  -  *4. If  Max  Taap  -  %4.1f\n", 
avg_cal_taaip,  tamp_nin>  tanp_max) ; 
braak;  “ 

caaa  ' 2' t 

poly_eolva() ; 
braak; 
caaa  '3' i 
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aet_conatant a ( } ; 
aet_atod_conatanta () ; 
break; 


> 


) 

}  while  (done 


caae  ' q' 
break; 

—  0); 


done  -  1; 


/*  poly.c  DTRC  Code  1543  J.  Blanton  5/89  */ 

/*  Program  to  determine  beat  valuea  of  alope  and  intercept  for  each 
probe  */ 

/•****•  Thia  file  containa  the  following  functional  ******* 

poly  aolve()  alope  and  intercept  for  3-coop  probea 

*/ 


♦include  "cal .h" 

PILE  *cal_flle; 

extern  double  hf_voltage_mean[30] [3] ; 
extern  float  avg_cal_teap; 
extern  double  pitot_velocity_mean(30] ; 
int  num_of_cal_pointa; 

poly_aolve ( ) 

( 

float  elope[3] , incep[3]  ,eoeff [3]  ; 

float  aum_hf [3] , aum_hf 2 [3] , aum_hf_pitot [3] ; 

float  aum_pitot,  aum_jpitot2;  ~ 

float  paeudo_vel [4] ; 

float  errl,err2,err3; 

double  aqr, temper; 

int  i,j,teat; 

float  num; 

char  check [81] ; 

double  top; 

float  n_inv,exp; 

get_cal_data ( ) ; 

num  -  num_of_cal_pointa; 

/*  LOOP  POP  DETERMXMXMB  BEST  1/M  */ 
printf("\n  exp  #1  #2  #3"); 

for <exp-l. ;exp<-2. /exp  +—  .05)  [ 
n_inv  -  1.  /  exp; 

/*  Zero  out  */ 
for (i— 0;  l<-2  ;  i44)  { 
aum_hf [i]  -  0. ; 
aum_hf2[i]  *0.; 

•urn  hf_pitot[i]  -  0 

) 

eum_j>itot  ■  0.;  aumjpitot2  —  0.; 

for  (i— 0;  i<-2  ;  i44)  ( 

for  (3-0;  3<num_of_cal_pointa;  j++)  { 

aumjhf  [£)  4-  (hf_voltage_mean  [  j )  [i  ]  * 
hf_yoltage_mean[ J) [i] ) ; 
au*_hf2 (i]  "4—  pow (hf_voltage_mean[ j] [i] , 4 . ) ; 
aum_hf _pitot[i]  4-  pow(hf_yoltage_mean[ j] [i] , 2. ) * 

pow (pitot  velocity  mean[jj,n  i nv); 

) 

» 

for  (J-0;  3 <nu»_o f _ca l_po i n t a ;  3++)  ( 

eum_pltot  4-  pow(pitot_velocity_mean[3J,n_inv); 
aum_pitot2  4-  pow (pitot  velocity  mean[3], (n  inv*2.)); 

] 
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/*  Find  Slope*  and  Intercepts  */ 

£or(i-0;  i<-2;  i++)  < 

sloped]  -  (sum_hf_pitot [i]  -  sura_pitot*sum_hf [i] /num) / 
(sum  h£2[i]  -  aum_hf[i] *aum_h£(i) /num) ; 

incap [1]  —  (aum_pitot-alopa[i]*sum_hf(lj)/num; 

top  -  sura  h£jpitot[i]-sum_h£[i;*sum_pitot/num; 

coaff [i]  -  top  *  top  /  " 

( (sum_hf2 [i]  -  sum_h£[i] *sum_hf (i) /num) * 
(sum_pitot2  -  sum_pitot*sum_pitot/num) ) ; 

printf ("\nexp-%5.3f  Cosff  -  %6.4f  %6.4f  %6.4f", 

exp, coaff [0] , coaff [1] , coaff [2] ) ; 


) 


/*  LOOP  FOR  DETERMININE  BEST  1/M  */ 
do  { 

print f ( " \n\nEntar  a  valua  of  n  (0  **  (1/n) !  ; 

scant ("%f", Sn_inv) ; 
n_inv  -I.  /  n_inv; 

/*  Zero  out  */ 
for (i«0;  l<-2  ;  i++)  ( 
sum  hf  [i]  —  0. ; 
sum_hf2 [i]  -  0.; 
sum_hf_pitot [ i ]  “  0 . ; 

) 

sum_pitot  -  0 . ;  sum_pitot2  -  0 . ; 

for (i-0;  i<-2  ;  i++)  ( 

for  (j-0;  j<num_of_cal_J>Oints;  j++)  ( 

aumhf [i]+-  <hf_voltage_mean[j]  dl* 
hf _voltage_mean [ J ] [1] ) ; 
sum_h£2 ti]  +-  p of (hf_vol t age_mean [ J ] [i] ,  4 . ) ; 
sum  hf  pitot [i]  +-  pow (h£_voltage_mean [ j ]  [i],2.)* 

“  pow  (pitot_valoclty_mean  [  j  ] ,  n_lnv) ; 

) 

> 

for  (j-0;  j<num_of  caljpoints;  j++)  ( 

sum_pitot  pow (pitot_velocity_*ean [ j] ,  n  lnv) ; 

sum_pitot2  +-  pow <pitot_velocity_mean [ j ] , (n_inv*2 . ) ) f 

) 

/*  Find  Slopes  and  Zntarcapts  */ 
printf ("\nn  -  *3 .3f\n", 1. /n_lnv) ; 

printf ("\n  Slops  Intarcapt  Corr  coaff\n"); 

for (i— 0;  i<— 2;  i++)  ( 

slops [i]  -  (sum  hf jpitot  [ij  -  Bum_pitot*sum_hfIi]/num)/ 
(sum_hf2[i]  -  sum_hf  [i]*sum_hf  (ij/num) ; 

incap [ i ]  -  (sum_pitot-slopa[i)*su*_hf[i])/num; 

top  -  sumjhf jpitot dj -sum_h£ [i J  •sum_pitot /num; 
coaffti]  -  top  *  top  / 

( (sum_hf2 [i]  -  sum_hf [i] *sum_hf (1) /num) * 
(sum_pitot2  -  sum_pitot*sum_pitot/num) ) ; 

printf ("\n#%d  %6.4f  %6.4f  *6.4f", 

i+1, slops (ij , inoap (i J , coaff [1] ) ; 


tast  -  If 

do  (  /*  start  of  plotting  do  */ 

printf ("\n\nHit  n  to  try  naw  value, 'p'  to  plot  or'q'  to  acaept  values  "); 

scanf  ("%s", check) ; 

switch  (tolower (check [0J)>  { 

'p's  /*  Plot  data  */ 

plot_eal  <hf_voltage_mean, 
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pltot_vslocity_ms«n, 

slop*, incap, n_inv, nuzn_of_cal_points) ; 
printf ("\nn  -  %5 . 3£\n", 1 . /n_inv) ; 

printf ("\n  slop*  —  Intercept  Corr  coeff\n") 

for  (i-0;i<-2;i++) 

printf ("\n*%d  %6.4f  %6.4f  %6.4f 

i+1, slope U] , inceptl] , coeff  [1] ) ; 

break; 
case  'q'  : 

test  *■  0  ; 

break; 
case  'n' ; 

teat  -  0; 
break; 

)  /*  End  of  switch  */ 

}  while  (test  —  1) ;  /*  End  of  plotting  do  */ 

}  while  (tolower  (check  [0] )  !•  'q'); 

printf ("\nDo  you  want  to  save  calibration  valuea  to  file?  (y/n)  ") ; 
scanf ("%s", check) ; 

if  (tolower  (cbeck[0] )  I-  'y')  { 

printf ("\nOne  more  chance-  Don't  you  really  want  to  aave?\n\ 

(Data  will  be  forever  lost)  (y/n)  "); 
scanf ("%e", check) ; 

) 

if  (tolower (check[0] )  —  'y')  { 

/*  Save  alopes,  intercepts,  and  n_inv  to  calibration  file  */ 
char  fil*_nams[37] ; 

sprintf  (file_name,  "/home/blanton/bf /ealdata/spdeal.  %s",probe_num)  ; 
cal  file-fopen(file  name,  "w") ; 

if  Tlcal_fila)  { 

printf ("***  %a  could  not  be  opened\n", file_name) ; 
exit  () ; 

) 

printf ("\n  Saving  calibration  file  %s  . . . . \n\n", file_naaa) ; 
fprintf (cal_f lie, "%5 . 3f  %5 . 2f \n», 1 . /n_lnv, avg_c*l_temp) ; 
for (1-0;  l<-2;  1++) 

fprintf (cal_file, »%6. 4f  %«. 4f\n", slope [1] , incep[i] ) ; 
f close (cal  file) ;  ~ 

) 

/*  Check  calibration  data  */ 

printf ("\n  Pseudo 41  Pseudo 4 2  Pseudo 4 3  Ref.  VelocltyNn") ; 

exp-1 . /n_lnv; 

for (i«0; l<num_of  cal_point s ;  1++ )  { 
for(j-0; j<3; J++)  ( 

sqr-hf_voltage_mean[i] [ j] ; 
tempor-incep [ jf  +  slope (j]*sqr*sqr; 
pseudo  vel[j]  -  pow (tenpor,  (double)  exp); 

) 

errl- (pseudo_vel [0] -pitot_veloclty_mean[i] ) /pitot_velocity_sman[i] *100. ; 
err2- (pseudo_vel [1] -pltot_velocity_mean [1] ) /pitot- velocltyjsean [i  j  *100 . ; 
err3- (pseudo  vel  [2]  -pitot  velocity  mean [1]) /pitot-' velocity- aman[ij*100. ; 
printf  ("#%d%10 . 4f  %11 . 4f%ll .  4f%lf.  4f*7 . 2f  %7 . 2f%f.  2f\n",  i+1, 
pseudo_vel[0] , pseudo_v*l [1 J ,pseudo_vel (2] , 
pitot  velocity  ms  an [T] ,errl,err2,err3) ; 

) 

) 

/*  data_in.c  */ 

/•*•*••  This  file  contains  the  following  functions;  ******* 

get_cal_data()  gets  calibration  data 


•include  "cal.b" 

TILE  *cal_vel; 

TILS  * constants, * output; 
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extern  int  num_of_cal_pointa ; 
extern  double  hf_voltage_mean[30]  [3] ; 
extern  double  pitot_veloclty_mean[30] ; 
extern  float  cal_temp[30] ; 
extern  char  probe_num[4] ; 
extern  float  avg_cal_temp; 
float  avg_cal_te«p_c; 
char  out_f ile_name [41] ; 

get_cal_data ( ) 

{ 

lnt  1; 

char  dummy [81] ; 

float  acratcb, acratchl,acratch2; 
float  teap_corr, cal_temp_c; 
double  aum; 
float  temp_max-0  ; 
float  tenp_min-100; 

/*  Find  number  of  calibration  polnta  */ 
char  £ ile_name [ 41 ] ; 

aprintf (file_name, "/hone/blanton/hf /caldata/rawapeed. %#" , probe_num) ; 
cal_vel-£open (flle_name,  "r"); 
if  (tcal_yel)  ( 

printfC***  speedcal.%a  could  not  be  opened\n",probe_num) ; 
exit(); 

) 

num_of_cal_polnt»«0 ; 
while  T I feof (cal_vel ) )  { 

fgets  (duagy,  80, cal_yel) ; 
num_of_cal_polnta++ ; 

> 

feloae (cal_vel) ; 
num_of_cal_polnte— ; 

/*  f a_velocity-dvector (1, num_of_cal_polnta) ; 
alpha-  vector  (1,  nuin_of_cal_point  a)  ; 
hf_yoltage-daatrlx  (T,  num_of_cal_pointa,  1,2); 
error-  vactor(l,num_of_cal_pointa); 
cal_teap-  vector (l,num_of_cal_pointa) ;  */ 

cal_vel-fopen (file_naam, "r") ; 

•ua»-0; 

for  (i»0;l<num_of  cal_pointa; 1++)  { 

f acanf (cal_vel, "%lf  «lf  %lf  *lf  %f», 

thf_voltage_meen[l] [0] , 
thf_voltage_mean[l] [1] , 

Chf_voltage~maan[ij [2] , 
cpltot  velocity  maan[i), 

*cal_tewp[l] ) ;  “ 

if  (cal  teap[l)  <  temp  min)  temp  win— cal  temp[i) ; 
if  (cal~teap[l]  >  temp_»ax)  tamp“max-c*l“temp[i] ; 
aum  4—  cal_tamp[l]; 

) 

feloae (cal_yel) ; 

/»  avg_cal_teap  -aum/num_of_cal_polnta;  */ 

/*  uae  firat  temperature  aa  cal  teaperatuxe  and  correct  aubaequent  voltagea  to  thia  cal 
temperature  */ 
avg_cal_t  aap  ■  cal_temp  [0] ; 
avg_cal_temp_o- (cal_teap[0]  -32 . )  /1 . 8; 
for  (T-0 ; i  <num  of- cal_polnta ; i++)  ( 

cal_teap_c  -  (cal_teap[l] -32. )/l.l; 

teap_corr  -  aqrt  ( (250 .  -avg_cal_teap_c)  /  (230 .  -cal_teap_c) )  j 
hf_voltage_mean[l] [0]  *-  teap_corr; 
hf_voltage~*ean[ij [1]  *-  teap_corr; 
hf~voltage~mean[ij [2]  •-  tamp  eorr; 

> 

print f ("\nhvg  temp-«4.1f  Min  Temp  -  *4. If  Max  Temp  -  %4.1f\n". 
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avg_cal _ tetnp,  tsrcp_mi n ,  temp_max) ; 


> 

/*  cal .  h  */ 

♦include  <stdio.h> 

< include  <mr.h> 

# include  "/dtrc/include/nrutil.h" 
♦  include  <Jaath .  h> 

♦define  PI  3.1415926 


/*  Local  parameters  (note  explicit  typing) :  */ 


♦define  HULL  0 
•define  8IZB0F0CA  1000 
♦define  HULLIHTRAY  -2 

lnt  nfrasms;  /* 

char  probe_num[4} ; 
float  ‘alpha, *phi; 
double  *pitot_volt; 
double  **hf_volt; 

/*  float  gain,offaet;  /* 

enua  (  AFF  ,  LTHT  )  wind_tunnel ; 

/• 


♦  of  aaa^les  per  channel  */ 

/*  Humber  for  specific  probe  */ 

anemometer  gain  and  offset  */ 
path  number  for  m  488  */ 


int  pathno; 
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APPENDIX  C 


THREE-COMPONENT  HOT-FILM  ANGLE  CALIBRATION  PROGRAMS 

(Note:  These  programs  are  not  complete  in  themselves;  only  functions  (subroutines) 
pertinent  to  hot-film  angle  calibration  are  included). 

/*  Program  angleo.c 

J.  Blanton  May  1989  Coda  1543  DTRC 

Baaad  haavily  on  Pat  Purtall's  AN0V0.7TN  angla  calibration 
program  May  1987  for  a  PDP  11/73 

To  estioute  tha  anglaa  tba  hot  film  alamanta  make  with  tha  coord, 
ayatam.  Trial  dlraction  coaina  data  ara  raad  and  itaratad  upon  to 
minimise  arror. 

Input  filaa: 

angval . (probe#)  angla  and  paaudo  valocitiaa 

dangla.dat  trial  dlraction  coainaa 

*/ 

(includa  "angcal.h" 

#lncluda  <atdlo.h> 

FILE  *cal; 
float  START  KAPPA; 
float  KND  KAPPA; 
float  KAPPA_INC; 

float  **dlr_cos;  /*  dlr_coa[3] [31  */ 

float  *actual_val;  /*  actual_val [3]  */ 

float  *v_veloclty;  /*  v_yelocity [3]  */ 

float  *coaputed_yel;  /*  cowputedjvelocity [3]  */ 
float  *q_squared;  /*  <j_«qu«rwd[3]  /  »/ 

float  **lnv_mat; 

float  **arr7  /*  computed  -  actual  velocity  [ang_cnt]  [3]*/ 

int  *kappa_dc, kappa_num_mln; 
anum  (  TBS  ,  HO  )  print_val ; 

maln() 

{ 

int  1, j,k,dc_num,kappa_cnt,kappa_num; 

int  deptbjBum; 

float  kappa_adn; 

float  *avg_err, *avg_aum_aqr; 

float  *kappa_var; 

float  get_angle_f lie ( ) , gat_trlal_num ( ) ;  /*  function  declarations  */ 

float  gat_trial_coa() ; 
char  check [ 81 ] ; 
char  flla_nama[41] ; 

float  *adn_aum_aqr; 

int  *mln  noa;  /*  dir.coa  #  whara  min  arror  occurs  */ 

char  duaaylSl]; 

print_val  -  MO;  /*  Disable  printing  of  velocity  values  */ 

/*  Determine  which  probe  to  use  •/ 

■et_constants  () ; 

/*  Read  in  Angle  velocity  data  with  alpha  and  phi  angles  */ 
get_angle_file () ; 

/•  Allocate  mmeory  and  unit-offset  the  matlx  for  dir.  cosines  */ 

dir_cos  matrix (1,3,1, 3) ; 

actual_vel-vector (1,3); 

v_velocity— vector (1,3); 

co«puted_vel -vector (1,3); 
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q_squared-vector (1,3); 
lnv_»at-iaatrix ( 1 ,  3 , 1 ,  3 )  ; 
err— matrix  (1,  3, 1,  ang_cnt)  ; 
avg_err«* vector (1,4); 
avg~aum_aqr-vector (1,4); 
min_num-i vector  (1,  4)  ; 
mi  n__aum_aqr— vector (1,4); 

depth_num- (RND_KAPPA  -  START_KAPPA) /KAPPA_INC; 
kappa~v»f“vvctor (1, 100) ; 
kappa_dc-l vector (1, 100) ; 

/*  kappa_dc-l vector (1,30) ;  */ 

/♦Read  In  number  of  direction  cosines  end  Initial  value  of  yaw  constant**/ 
get_trial_num()  ; 

/♦  Loop  through  several  kappas  */ 
kappa_num-0; 

for (kappa-START_KAPPA;kappa<-END_KAPPA;kappa+-  XAPPA_IHC)  { 
kappa_num++ ;  _ 

for (i-1; i<-4 ; 1++) 

min_sum_sqr[i 1-1000; 

/*  Loop  through  calculation  for  the  number  of  direction  cosines  */ 
for (i-1; i<-dc_cnt ; i++)  ( 

/*  Read  in  a  set  of  trial  direction  cosines  »/ 
get  trial  cos (dir_coa, 1) ; 

/*  ~ 

printf ("\n%f  *f  %f  %f  %f  %f  %f  %f  %f\n", 

dir_cos[l] [1] ,dir_eos[ll [2J ,dir_cos[l] [3J , 
dir  coa[2] [lJ,dir_cos[2) [2] , dlr_cos[2] [3], 
dir_cos [3] [1] ,  dir  cos[3] [2], dir  cos[3][3]); 

*/ 

compute_velocity (i) ; 

/*  Compute  statistics  */ 
for(k-l;k<-3;k++)  ( 

sum_square (err [k] , ang_cnt, fiavg_sum_aqr [k] ) ; 

/*  stats (err [k] , ang_cnt, «avg_err [kj , savg_sum_sqr [k] ) ;  */ 

if (avg_sum_sqr[k]  <  min_sum_sqr [k] )  ( 

~min_sum_sqr[kl  avg^sum_sqr  [k]  ; 
min_num[k)  -  i; 

1 

) 

/*  Compute  average  of  averages  */ 

avg_err[4)- (avg_err [1 J  +  avg_err[2]  +  avg_err[3] ) /3. ; 
avg_sum_eqr[4]-(avg_sum_aqr[ll  + 

avg_sum_aqr[2]  + 
avg_sum_sqr [3 j )/3 . ; 
if  (avg_sua_sqr[4]  <  min_sum_sqr [4] )  ( 

min_sum_sqr[4]  -  avg_sum_sqr [4] ; 
min_num[4]  -  i; 

1 

/•printf ("\nCos#%d  Ovar-%5.2f  War-%5 . 2f  Wvar— %5 . 2f  Avg  var— %5.2f" 

,1, avg_sum_sqr [1] , avg_eum_eqr [2J , avg_aum_sqr[3J , avg_sum_sqr [41 ) ;  */ 

1  /*  Xnd  of  dir.  cos.  loop  */ 

/*  print  out  minimum  errors  and  coresponding  direction  cosines  */ 
printf ("\nMaxinum  alpha  -  %4.1f  Xappa  -  %6.4f\n", 
alpha [ang_cnt] , kappa) ; 

/*  for(k-l;k<-3;k++)  (  */ 

/*  printf ("\nMin  avg  var  error#%d  -%5.3f  dir.  cos.  #%d" 

,k,min_sum_sqr[kj ,mln_num[k) ) »  */ 

/*  get_trlal  cos <dir_cos,sdn_num[k] ) ;  */ 

/*  print_dc(3ir_cos)7  */ 
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/*  >  */ 

printf ("AnMin  avg  var  •rror  —  %5.3f  dir.  cos.  #%d\n\n" 

, min_sum_sqr [ 4 ] , mln_num [ 4 ] ) ; 
gst_trial_cos (dir_coa, min_num[4] ) ; 
print_dc (dir_eoa) ; 

kappavar [kappa_num] -min_aum_aqr [ 4] ; 
kappa_dc [kappa_nua] •min_num [ 4 ] ; 

)  /*  End  of  kappa  loop  */ 

kappa_num_min-l  ; 
kappa_min-8TART_KAPP A ; 
kappa_num— 0 ; 

for <kappa-START_KAPPA; kappa<-BND_KAPPA; kappa+-  KAPPA_IHC)  ( 
kappa_num++; 

printf ("\nKappa  -%5.2f  Min  tnaan  square  arror  -  %5.3f  dir.  cos.  #%d" 
,  kappa, kappa_var [kappa_num] , kappa_dc [kappa_num] ) ; 
if (kappa_var [kappa_num]  <  kappa_var [kappa_num_min] )  { 
kappa_num_mi n-kappa_num ; 
kappa  min-kappa; 

> 


) 

printf ("\n\nKappa  »%5.2f  Min  Min  naan  aquara  arror  -  %5.3f  dir.  coa.  #%d" 
, kappa_mln, kappa_var [kappa_num_nin] , 
kappa_dc[kappa_num_min] ) ; 

gat _t rial  cos (dir_coa,kappa_de[kappa_num  min]); 
printf ("\n"); 
prlnt_dc (dir_coa) ; 
kappa— kappa_min; 

printf ("\nDo  you  want  to  aava  calibration  valuaa  to  fila?  (y/n)  ") ; 
acanf ("%a", chack) ; 

if  (tolowar  (chack[0] )  'y' )  { 

printf ("\nOna  mora  chanca-  Don't  you  raally  want  to  aava?\n\ 

(Data  will  ba  foravar  lost)  (y/n)  "); 
acanf ("%b", check) ; 

) 

if  (tolowar  (chack [0] )  —  'y')  { 

sprint f (fila_naaa, "/boma/blanton/hf /caldata/dircoa. %a",probe_num) ; 
printf ("\n%e?, fila_nama) ; 
cal-fopen(file_nane,  "w"); 
if  (leal)  ( 

printf ("***  %a  could  not  ba  opanad\n", file_naae) ; 
exit  () ; 

) 

f printf (cal, "%4 . 2f \n", kappa_min) j 

get_trial_coa (dir_coa, kappa_dc[kappa_nun_nln] ) ; 

for  (i-l;i<-3;i++)  { 

for ( J»l; j<-3; j++)  ( 

fprintf (cal, "%6. 4f  »,dir  coa[iJ[jJ>; 

)  fprintf (cal, "\n"); 

) 

fcloaa (cal) ; 

) 

/*  Chack  calibration  data  */ 
chack  cal  () ; 


compute_velocity (i ) 
int  i; 

( 

int  j,k; 

if  (print_wal  —  TES  )  ( 

printf ("\nkappa  -  %4. 3f", kappa) ; 

printf ("\n  al  phi  ral  act  coa?  %%arr  *); 
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} 

/*  Loop  through  alpha  and  phi  angle  val  data  computing  arrora  */ 
for ( j-1 ; j<-ang_cnt ; j++)  { 

/*  Compute  U,v,N  (actual  velocity  in  probe  coordinatea) 
given  alpha  and  phi  */ 

actual_vel[l]-ref_vell j]*coa (alpha [ j] *PI/180 . ) ; 
actual_vel 12] »ref_vel [ j] *aln (alpha [ j) *PI/180 . ) * 

coa (phi I  j] *PI/180 . ) ; 

actual  vel[3]-  -1 . *ref_vel [j] *ein (alpha [j] *PI/180 .) * 

ain(phi[jl*PI/180. ) ; 

/*  Compute  Q' a  given  assumed  direction  coainea 
and  calibration  velocity  */ 

/*  if  (print_vel  —  YES  ) { 

print f  ("\nkappa-%4 . 3 f”, kappa)  ; 
printf ("\ndir  coaU) [i] —  %6.5f  %6.5f 
%6. 5f ", dir_coa [1J [1] ,  dir  coa[l) (2) , dir  co«[l)[3]); 

)  */ 

for(k-l;k<-3;k++)  { 

q_squared[k] -pseudo ( j] (k) *pseudo [ j) [k]* 

(1.  +  (kappa*kappa  -  1.)* 
dir  coa tl J Ik] *dir_coa II] (k) ) ; 

) 

/*  Compute  Ef fectiveneaa  Matrix  (Kappa  aquared)  and  invert  */ 
kappa_matrix(inv_mat) ; 

/*  Compute  VJ  by  multipling  Q*2  by  inverted 
kappaA2  matrix  */ 

for (k«l ; k<“3 ; k++)  ( 

v_velocityIk]-inv_mat(k]  (1) *q_aquared[l]  + 

inv_mat (k] [2] *q_aquared(2]  + 
inv_mat [k] [3]  *q_aquared[3] ; 
if  (v_velocity[k]  <  0.)  ( 

printf ("**  Negative  V  **  "); 
printf ("alpha-%4. If  phi-%5.1f  i-»d\n", 
alpha  [  j  ] ,  phi  t  j  I » i ) ; 
v_veloclty [k]-  -l*v_velocity [k] ; 

} 

v_velocity Ik] -aqrt (v_velocity Ik] ) ; 

) 

/*  Compute  u,  v,w  (computed  velocity  in  probe  coordinatea) 
given  direction  coainea  and  Vj' a  */ 

for (k—1 ; k<“3 ; k++ ) 

c  ompu t  ed_ve 1 Ik] -dir_coa Ik] 11] *v_yelocity II]  + 

dir_coo Ik] 12] *v_ve locity {2]  + 
dir_coa[k] I3]*v_velocltyI3] ; 

/*  Compute  error  (difference  bewteen  actual 
and  computed  velocity)  */ 
for(k-l;k<-3;k++)  { 

•rr Ik] [ j ] -computed  vel Ik] -actual  vellk]; 

) 

if  (print  vel  — —  YES  )  ( 

printf ( -\n%5 . If %6 . If %6 . 2f  %7 . 2f%7 . 2f %5 . If %7 . 2f %7 . 2f%5 . If %7 . 2f %7 . 2f %5 . If " , 

alpha [ j ] , phi I j ] , ref_vel l j ] , 

actual_vel (1] , computed_vel [1 J , err II] I J] /actual_vel II ] *100 . , 
actual~vel(2] ,computed~vel[2] ,err[2] [ j]/actual~vel[l]*100. , 
actual_vel (3) , computed  vel 13] , err [3] [ j ] /actual- vel [1] *100 . ) ; 

} 

)  /*  end  of  alpha-phi  loop  */ 

if  (print_vel  —  YES  ) 
printf ("\n"); 


check  cal  <) 

( 

get_trial_coa (dir_coa, kappa_dcfkappa_num_min] ) / 
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print_val  -  YES; 
computs_vslocity (0) ; 


print_dc (dlr_co«) 
float  **dir  coa; 

< 

int  i,J; 

float  Bum; 

for  (i-l;i<-3;l++)  { 

for(j-l; j<-3; j44)  { 

prlntf  ("%6.  4f  ",  dir_eosflH  j)  > ; 
)  prlntf <"\n"); 

) 

prlntf ("\nSumA 2  Col."); 

for(i-l;i<-3;144)  { 

sun-  dlr_coa[l] [l]*dlr_coa[lJ tl]  + 
dlr_cos[2]  [l]*dlr  coa[2]tl]  + 
dir_coa[3) [1] *dir~cos [3] til ; 
prlntf <"  #%d  -  %f",i,sum); 

) 


prlntf ("\nSumA2  Row  "); 
for(i-l;i<-3;l++)  { 

bub-  dlr_coa [1] [1] *dlr  cos[l] tl]  + 

dir_cos[l] [2)*dir-eos[i] [2]  + 
dir_coa [i] [3]*dlr“cos[U (3); 
prlntf ("  #%d  -  %f",i,sumj; 

) 

prlntf ("\n") ; 


> 

sum_squara (x, nun, sum_sqr ) 
float  *x,  *sum_aqr; 
lnt  nun; 

( 

lnt  1; 

doubla  sum-0 ; 
f or (1-1; l<-num; 14+)  { 
aum  4-  x[l]*x[l] ; 

) 

*sum_sqr  -  (float)  sqrt(sum); 

> 

/*  data_in.c  DTRC  Coda  1543  5/89  */ 

/•*•**•  This  tlla  contalna  tha  following  functions:  ******* 

gat_angla_fila()  gats  angla  calibration  data 

gat~trial_num()  gats  #  of  diraction  coslnaa 

and  kappa  and  kaapa  f 11a  opanad 
to  raad  In  valuaa  as  naadad 

gat_trial_coB ()  Assign  trial  diraction  cosinas 

~  from  tamp  array  */ 


finduda  "angcal.b" 

/**•**  Raad  In  Angla  valodty  data  with  alpha  and  phi  anglas  */ 
gwt_angl*_flla  () 

( 

lnt  1; 

char  dummy (81); 

/*  rind  numbs r  of  angla  calibration  points  */ 
char  fil*_na»s (37] ; 

sprintf (fila_naaw, "/homa/blanton/hf/caldata/angwal. »s",proba_num) ; 
ang_wal-fopan (flla_nsms,  "r "); 

If  ( I ang_T*i )  ( 

prlntf  ("***  AMOVE!.  %s  could  not  ba  opanad\n" ,  proba  num) ; 
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•xit ( ) ; 

> 

ang_cnt-0 ; 

while  < I feof (ang_yel) )  { 

fgets (durany, 80, ang_vel) ; 
ang_cnt++; 

) 

fclose (ing  v«l) ; 
ang_cnt — ; 

/*  unit-offset  alpha  and  phi  (aaa  Numerical  Recipes  in  C  p.l5sl6)  */ 
alpha-  vector (1 , ang_cnt ) ; 
phi-  vector (1, ang_cnt) ; 
ref_vel-  vector (1 , ang_cnt ) ; 

pseudo-matrix (1, ang_cnt, 1, 3) ;  /*  Numerical  Recipes  matrix  allocation  */ 

ang_vel— f open (f ile_name, "r") ; 
for  (i-l;i<-ang_cnt;i++)  { 

fscanf <ang_vel, "%f  %f  %f  *f  %f  %f ", ialpbs (1) , tphi [i] , tpseudo [i] [1] , 
t pseudo [i] [2] ,  /*  switched  2  and  3  */ 

tpseudo [i] [3] , tref_vel [i] ) ; 

/*  print f ("\n*f  %f  *f  %f  %f  %f",alpha[i] ,phi[i] ,pseudo[i) [1] , 
pseudo [i] [2] , pseudo [i] [3] , ref_vel [i] ) ;  •/ 

) 

fclose (ang_vel) ; 

} 

/*  Read  in  number  of  trial  direction  cosines  and  yaw  constant  */ 

/*  Keep  file  open  to  read  in  values  as  needed  */ 

get_trial_num ( ) 

< 

int  i , j , dc_num; 

d_angle-f open ("dangls.dat",  "r") ; 
if  (ld_angle)  { 

printf ("***  DANOLS.DAT  could  not  be  opened") ; 

•xit(); 

) 

fscanf (d_angle, "%d  %f ", £dc_cnt , tkappa) ; 
printf  (*4d  *f \n", dc_cnt, kappa) ; 
dc_terap-matrix  (1 ,  dc_cnt  ,1,9); 
for(i-l;i<-dc_cnt;i++)  { 

for (j-l; J<-9; J-M-)  { 

fscanf (d  angle, "%f", sdc  temp[i][jl); 

) 

fscanf (d  angle, "%d" , <dc_num) ; 

) 

fclose (d_angle) ; 

> 

get_trial_cos (dc, nun) 
float  **dc; 
int  num; 

( 

/*  Assign  trial  direction  cosines  from  teqp  array  */ 

int  j,k,dc_num,  ent; 
ent— 1; 

for ( J-l; J<— 3; j++)  [ 

for (k— l;k<-3;k++)  { 

dc[j][k]-dc  temp  [num]  [cnt++]; 

) 

) 

) 

/*  movagv . c  DTRC  Coda  1343  5/89  V 

/*  This  file  contains  : 

state (v,n,avg,var)  Coapute  the  mean  and  variance  of  a  vector 

of  data  (v) 
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where : 


v  data  array  (must  be  unit -off aat) 

n  n unbar  of  points 

avg  Man  value  returned  (must  be  pointer) 

var  variance  value  returned  (must  be  pointer)  */ 


♦include  "nrutll.b" 

stats  (v,  n,  avg,  var) 
float  *v,*avg,*var; 
int  n; 

( 

int  i; 

float  avgold; 
for(i-l;i<-n;i++)  ( 
if(i— 1)  { 

*avg— v[l J ; 

* var— 0 . j 
avgold-  *avg; 

) 

if  (1—2)  ( 

* avg- (avgold  +  v[2])/2.; 

*var-0. 5* (avgold  -  v[2] )* (avgold  -  v[2J); 
avgold-  *avg; 

) 

if(ll-l  C  il-2)  { 

*avg  -  ( (i-1) *avgold+v [il ) /i; 

*var  -  ((i-2)*  (*var)  +  (i-1) * avgold* avgold 
-  i*(*avg)*(*avg)  + 
v[i]*v(i])  /  (i-1.); 

avgold  «  *avg; 

} 

) 

) 

/*  angcal.b  */ 

♦include  <atdio.b> 

•include  "/dtrc/lndude/nrutil -h" 

♦include  <math.h> 

♦define  PI  3.1415926 
FILE  »ang_vel, *d_angle; 

float  kappa; 
int  dc_cnt; 
int  ang_cnt; 

float  ‘alpha, *phl, *ref_vel; 
float  **pseudo; 
float  **dc_temp; 
char  proba_num[4) ; 


/*  ♦  of  direction  cos-counter  */ 
/*  ♦  of  angle  cal  pts-counter  */ 

/*  alpha [ 1 , ang_cnt ]  etc  */ 

/*  pseudo [l,ang_cnt] [1,3)  */ 

/*  dc_te*p[l,dc_cnt)  [1,9)  */ 
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(This  page  intentionally  left  blank) 
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